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ABSTRACT 

We present the rest-frame colors and luminosities of ^ 25000 m« < 24 galaxies in the redshift range 0.2 < z < 
1.1, drawn from 0.78 square degrees of the COMBO-17 survey. We find that the rest-frame color distribution of 
these galaxies is bimodal at all redshifts out to z ~ 1 . This bimodality permits a model-independent definition 
of red, early-type galaxies and blue, late-type galaxies at any given redshift. The colors of the blue peak become 
redder towards the present day, and the number density of blue luminous galaxies has dropped strongly since z ~ I . 
Focusing on the red galaxies, we find that they populate a color-magnitude relation. Such red sequences have been 
identified in galaxy cluster environments, but our data show that such a sequence exists over this redshift range 
even when averaging over all environments. The mean color of the red galaxy sequence evolves with redshift in 
a way that is consistent with the aging of an ancient stellar population. The rest-frame Z?-band luminosity density 
in red galaxies evolves only mildly with redshift in a A-dominated cold dark matter universe. Accounting for the 
change in stellar mass-to-light ratio implied by the redshift evolution in red galaxy colors, the COMBO-17 data 
indicate an increase in stellar mass on the red sequence by a factor of two since z ~ 1 . The largest source of 
uncertainty is large-scale structure, implying that considerably larger surveys are necessary to further refine this 
result. We explore mechanisms that may drive this evolution in the red galaxy population, finding that both galaxy 
merging and truncation of star formation in some fraction of the blue, star-forming population are required to fully 
explain the properties of these galaxies. 

Subject headings: galaxies: evolution — galaxies: general — galaxies: luminosity function — galaxies: 
elliptical and lenticular — galaxies: stellar content — surveys 



I. INTRODUCTION 

The evolution of early-type galaxies since early in the Uni- 
verse's history is a highly active area of research. The cur- 
rent world model for the evolution of the Universe is the Cold 
Dark Matter paradigm; the incarnation which appears consis- 
tent with most observations adopts ft^ = 0.3, JIa = 0.7, and 
Ho = 100/ikms"' Mpc"', where h^Q.l (the ACDM paradigm; 
e.g., Freedman et al. 2001; Efstathiou et al. 2002; Pryke et al. 
2002; Spergel et al. 2003). In this model, early-type galax- 
ies form hierarchically through mergers of pre-existing galaxies 
(e.g.. White & Frenk 1991; Barnes 1992; Cole et al. 2000). A 
generic prediction of this type of model is an increase in the 
stellar mass density of the early-type galaxy population since 
1 , as large galaxies are assembled primarily at relatively late 
times (e.g., Aragon-Salamanca, Baugh, & Kauffmann 1998; 
Kauffmann & Chariot 1998). Therefore, an important test of 
the hierarchical formation scenario is to quantify the evolution 
of early-type galaxies since z ~ 1 . In this paper, we discuss 
the color distributions of a sample of ~ 25000 galaxies from 
the COMBO-17 survey ('Classifying Objects by Medium-Band 
Observations in 17 Filters'; Wolf et al. 2003), focusing on an 
objective definition of early-type galaxies at all redshifts. We 
then use this objectively-defined sample of early-type galaxies 
with 0.2 < z < 1.1 to examine the evidence for a substantial 
build-up of stellar mass in red galaxies since z ~ 1 . 



In the local Universe, early-type galaxies have historically 
been defined in terms of morphology. Loosely speaking, galax- 
ies with dominant smooth spheroidal components, and at most 
a modest and smooth disk, will be classified as an 'early-type' 
galaxy. These galaxies tend to be red in optical color both lo- 
cally (e.g., Schweizer & Seitzer 1992) and out to z ^ 1 (e.g., 
Kodama, Bower, & Bell 1999). This definition of early-type 
galaxies is being challenged by the recent demonstration that 
galaxies are distributed in color space in a bimodal distribution 
(see, e.g., Strateva et al. 2001; Hogg et al. 2002; Blanton et 
al. 2003a). One peak is red, and consists mostly of non-star- 
forming galaxies earlier than Sa in morphological type. The 
other peak is blue and consists primarily of star-forming galax- 
ies later than Sb in morphological type (Strateva et al. 2001). 
Because of this bimodality in color and therefore star forma- 
tion history (SFH), it is defensible, and perhaps more natural, to 
define early types in terms of colors (i.e., spectral energy distri- 
bution or SFH). These different definitions seem largely consis- 
tent (e.g., Schweizer & Seitzer 1992; Strateva et al. 2001; Bell 
et al. 2003, 2004). In this paper, we split the galaxy population 
into red and blue subsets based on their rest-frame optical col- 
ors, exploring the evolution of these populations with cosmic 
time. 

In the local Universe, the blue-peak galaxies (later types) 
show a scattered but systematic variation of color with mag- 
nitude, in the sense that luminous galaxies tend to be some- 



2 



Bell et al. 



what less blue (e.g., Tully et al. 1998; Hogg et al. 2002). These 
color changes, when examined in detail, are due to changes 
in the mean ages, metallicities, and dust contents of galax- 
ies with galaxy magnitude, such that brighter galaxies tend to 
be older, dustier, and more metal-rich (see, e.g., Tully et al. 
1998; Bell & de Jong 2000). In contrast, red-peak galaxies 
(earlier types) have a tighter and well-defined relationship be- 
tween color and magnitude, such that bright galaxies are typi- 
cally redder (e.g., Sandage & Visvanathan 1978; Bower, Lucey, 
& Ellis 1992; Schweizer & Seitzer 1992; Terlevich, Caldwell, 
& Bower 2001). This color-magnitude relation (CMR) is well- 
estabUshed in overdense and cluster environments, where early- 
type galaxies are much more common (Dressier 1980; Dressier 
et al. 1997). Nevertheless, the CMR also exists among present- 
day field early-type galaxies (e.g., Sandage & Visvanathan 
1978; Schweizer & Seitzer 1992; Hogg et al. 2003). The scat- 
ter, slope and evolution of the CMR in clusters of galaxies are 
all consistent with the interpretation of the CMR as primar- 
ily a metallicity sequence of old galaxies, where more massive 
galaxies are more metal-rich (e.g.. Bower, Lucey, & Ellis 1992; 
Kodama & Arimoto 1997; Vazdekis et al. 2001; Bernardi et al. 
2003). A modest age spread is possible (Trager et al. 2000). 

The evolution of the properties of individual early-type galax- 
ies with redshift has been studied in clusters of galaxies, where 
the assembly of sizeable galaxy samples has been feasible. It 
appears as if individual early-type galaxies were formed rela- 
tively quickly at high redshift and simply aged to the present 
day through so-called passive evolution (e.g., Kodama & Ari- 
moto 1997; Kelson et al. 2001; van Dokkum et al. 2003). This 
has been used to justify 'monoUthic collapse' models of early- 
type galaxy formation, where these galaxies form quickly early 
in the history of the Universe, and simply age to the present 
day (e.g., Eggen, Lynden-Bell, & Sandage 1962; Larson 1975; 
Kodama & Arimoto 1997; Nulsen & Fabian 1997). Yet, mono- 
lithic collapse, at best, is mildly unphysical; the hierarchical 
build-up of galaxies through mergers and interactions must hap- 
pen at some level (e.g. White & Frenk 1991). 

An elegant interpretation of the apparent contradiction be- 
tween seemingly passive evolution and the hierarchical build- 
up of early-type galaxies has been put forward by van Dokkum 
& Franx (2001, progenitor bias). They note that studies of the 
properties of early-type galaxies in isolation (rightfully) fail to 
study the star-forming precursors of galaxies that are later des- 
tined to become early types. Thus, individual early-type galax- 
ies appear old at all redshifts regardless of whether they were 
formed in a burst at redshift infinity or whether their forma- 
tion happens through violent mergers which consume all of the 
gas in the interacting galaxies, leaving a spheroidal remnant 
(e.g., Toomre & Toomre 1972; Barnes & Hernquist 1996). This 
is a common theme in many contemporary studies of early- 
type galaxies; the distinction between when the stars in early- 
type galaxies were formed, and when the galaxies were assem- 
bled into recognizable entities. Studies of galaxies selected by 
strong gravitational lensing largely avoid this bias by selecting 
by mass (e.g., Kochanek et al. 2000; Rusin et al. 2003; van de 
Ven, van Dokkum, & Franx 2003); yet, these studies also find 
rather ancient-looking stellar populations. 

Thus, a more important constraint on the evolution of early- 
type galaxies — and therefore, a stronger test of galaxy for- 
mation models — is the evolution of their number density and 
luminosity function (whether the early-type galaxy population 
is defined by morphology or color). Here, the situation is some- 



what less clear Many studies find that the evolution of the 
early-type population is close to passive up to z 1, with per- 
haps a slight decrease in number density above z ~ 1 (e.g., Lilly 
et al. 1995; Schade et al. 1999; Lin et al. 1999; Cimatti et al. 
2002; Firth et al. 2002; Im et al. 2002; Chen et al. 2003; Pozzetti 
et al. 2003). Other studies find relatively rapid evolution. For 
example, Kauffmann et al. (1996) find number evolution of a 
factor of three by z ^ 1, although Totani & Yoshii (1998) argue 
that spectroscopic incompleteness causes this apparent evolu- 
tion, and state that Kauffmann et al's results are consistent with 
passive evolution. Wolf et al. (2003) find a rapid increase in the 
luminosity density of galaxies with the colors of present-day 
early-type galaxies from z ^ 1 until the present day; yet, pas- 
sive evolution and reddening of the galaxy population is likely 
responsible for much of this. 

Nevertheless, a problem shared by most surveys of early- 
type galaxies at z > 0.5 has been susceptibility to large scale 
structure. Surveys to date have typically had samples of only 
50-1500 early-type galaxies (e.g., Kauffmann & Chariot 1998; 
Schade et al. 1999; Lin et al. 1999; Firth et al. 2002; Chen et 
al. 2003; Pozzetti et al. 2003). Furthermore, these samples are 
usually from relatively small areas on the sky (the largest are 
CN0C2 and the LCIRS with 1400 square arcminutes each; Lin 
et al. 1999; Chen et al. 2003). Worse still, early-type galax- 
ies are strongly clustered (e.g., Dressier 1980; Dressier et al. 
1997; Moustakas & Somerville 2002; Daddi et al. 2002), which 
means that the measured number evolution is particularly sen- 
sitive to large-scale structure. Thus, to accurately measure the 
evolution of the stellar mass in early-type galaxies, it is neces- 
sary to use large samples drawn from large areas on the sky.' 

In this paper, we bring the largest present-day sample of 
galaxies with 0.2 < z < 11 to bear on this important prob- 
lem. The COMBO- 17 survey has to date imaged and fully ana- 
lyzed three ~0.25 square degree fields in 5 broad-band and 12 
narrow-band filters (data for a fourth field has been taken but 
is not yet analyzed). The 5 broad and 12 narrow-bands allow 
the construction of photometric redshifts and rest-frame col- 
ors with extraordinary precision (one can think of the data as 
a /? ~ 10 ultra-low resolution spectrum rather as a photomet- 
ric dataset), yielding galaxy redshifts accurate to a few percent 
for a total sample of ^ 25000 galaxies in 0.78 square degrees 
(2800 square arcminutes) with niR < 24. We use this sample 
to discuss the colors of galaxies with 0.2 < z < 1-1, focusing 
particularly on the reddest galaxies at all redshifts. We put for- 
ward an objective definition of early-type galaxies, and use this 
definition to explore the evolution of the red galaxy population 
from redshift unity to the present day. 

This paper is set out as follows. The COMBO- 17 survey and 
the main sources of error are briefly discussed in §2. In §3, 
we explore the colors of the entire COMBO- 17 galaxy popula- 
tion with 0.2 < z < 1.1. Then, in §4, we focus on the reddest 
galaxies at all redshifts, exploring their CMR in more detail. 
Defining early-type galaxies as galaxies on the CMR, we then 
go on in §5 to explore the evolution of the luminosity function 
of red galaxies. We discuss the results in §6, and present the 
conclusions in §7. Throughout, we assume J7m = 0.3, f^A = 0.7, 
and H() = 100/zkms"' Mpc"'. In order to estimate the passive 

' The cosmic variance error decreases by a factor of two for more tlian an 
order of magnitude increase in comoving volume probed (if the area is contigu- 
ous), or for a factor of four increase in the number of independent pointings 
(Somerville et al. 2004). Therefore, to derive an accurate cosmic average a sur- 
vey with many independent pointings is ideal. Typically, surveys cover a few 
independent pointings for practicality; ours is no exception. 
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evolution of stellar populations as a function of lookback time, 
we adopt h = 0.7 to be consistent with the Hubble Space Tele- 
scope Key Project distance scale (Freedman et al. 2001) and the 
recent results from the Wilkinson Microwave Anisotropy Probe 
(Spergel et al. 2003). 

2. THE DATA 

The COMBO- 17 data, the sample selection, redshift estima- 
tion, construction of rest-frame luminosities, and completeness 
are all discussed in detail by Wolf et al. (2003, W03 hereafter). 
Wolf et al. (2001a, WOl hereafter). Wolf et al. (2001b), and 
in a forthcoming technical paper (Wolf et al., in preparation). 
Here, we briefly discuss the relevant hmitations and sources of 
uncertainty in the COMBO-17 data. 

2.1. Data and Sample Selection 

To date, the COMBO-17 project has surveyed 0.78 square 
degrees to deep limits in three pointings of the Wide Field 
Imager (WFl; Baade et al. 1998, 1999) at the Max Planck 
Gesellschaft/European Southern Observatory 2.2-m telescope 
at La Silla, Chile. COMBO-17 targeted the Chandra Deep 
Field South, an equatorial field, and a field centered on the 
z = 0.16 AbeU 901 cluster (for the coordinates of the three 
fields, see W03). We study the galaxy population at z > 0.2, 
therefore the presence of the low-z Abell 901 cluster of galax- 
ies in one of the survey fields will not significantly bias our 
analysis. The WFl has eight 2kx4k CCDs, a field of view of 
34' X 33', and apixel scale of 0."238/pixel. A total of- 160ksec 
per field were taken in 5 wide and 12 medium passbands with 
3640A< Aobs < 9140A; an ultra-deep 20ksec R-band with see- 
ing below 0."8 and a IOct limiting magnitude of 25.2, and deep 
UBVI and 12 medium passbands designed to allow the con- 
struction of low-resolution R ^ 10 spectra, or equivalently ac- 
curate photometric redshift estimates, for a total of ~ 25000 
objects with 0.2 < z< 1-2. The comoving survey volume is 
^ 10^/z"^Mpc^ in the interval 0.2 < z < 1.1, split between three 
disjoint roughly equally-sized fields. 

Galaxies were detected in the deep 7? -band frames using SEx- 
tractor (Bertin & Arnouts 1996), and we adopt the MAG-BEST 
magnitude as our total magnitude estimate (a Kron 1980, mag- 
nitude in uncrowded regions). Bertin & Arnouts (1996) show 
that this magnitude underestimates the flux in galaxies by ~6%. 
We do not correct for this effect as almost all local galaxy 
surveys adopt similarly biased Kron or Petrosian magnitudes 
(e.g., Blanton et al. 2003b; Skrutskie et al. 1997). The spec- 
tral shapes for /?-band detected objects were measured by per- 
forming seeing-adaptive, weighted aperture photometry in all 
17 frames at the position of the /?-band detected object using 
the package mpiaphot (Roser & Meisenheimer 1991, Meisen- 
heimer et al., in preparation). This package measures the peak 
surface brightness of all of the images smoothed to identical 
seeing (~ l."5) to maximize signal-to-noise, allowing one to 
construct the best possible photometric redshifts. Thus, the col- 
ors are essentially total for distant objects and are central colors 
for more extended, nearby galaxies. The effects of this aper- 
ture bias are mild, and are discussed later in §4. Photometric 
cahbration was achieved using spectrophotometric standards in 
each COMBO-17 field. All magnitudes are quoted in Vega- 
normalized magnitudes. The median galaxy in our sample with 
niR ~ 22 has observed galaxy magnitude errors of Sr — 0.01 
mag, 6b,vj ~ 0.05 mag, Su ~ 0.3 mag, and (^Medium for the 
medium passbands increasing from 0.03 mag at the red end to 
0.2 mag at the blue end. 



2.2. Photometric Redshifts, Spectral Templates and 
Rest-Frame Colors 

The full survey dataset is photometrically classified into stars, 
galaxies and AGN using the 17-passband 'fuzzy' spectrum 
(WOl). No attempt is made to use morphological information, 
as, for example, double stars may contaminate the galaxy cata- 
log, and compact galaxies may contaminate the stellar or AGN 
catalogs. For galaxy classification, we use Pegase (see Fioc 
& Rocca-Volmerange 1997, for an earher version of the model) 
model spectra. The template spectra are a two-dimensional 
age/reddening sequence, where a fixed exponential star forma- 
tion timescale t = 1 Gyr is assumed, ages vary between 50 Myr 
and 10 Gyr and the reddening E(B-V) can be as large as 0.5 
mag, adopting a Small Magellanic Cloud Bar extinction curve^. 
These spectra are defined in the interval 1216A< Aobs < 3^m 
for elliptical through to starburst galaxies. Trustworthy tem- 
plates shortwards of 1 2 1 6 A have not yet been included, Umiting 
redshift measurements at this stage to z < 1.40. For our subse- 
quent analysis the redshift range is further limited to z < 1 . 1 in 
order to have at least two filters redwards of the 4000A break. 
The template Ubrary does not evolve, and cannot perfectly de- 
scribe spectral types which are highly unusual in the local Uni- 
verse, such as post-starburst E+A galaxies. Furthermore, low- 
luminosity or type 2 AGN will be misclassified as galaxies, 
sometimes leading to increased redshift errors. For galaxies 
and AGN, galaxy redshifts are estimated simultaneously with 
the template type (WOl). A description of an earlier version of 
this methodology was presented by W03, and a full description 
and analysis of these photometric redshifts and classifications 
wiU be presented by Wolf et al. (in preparation). 

The galaxy redshift estimate quality has been tested by com- 
parison with spectroscopic redshifts for many hundreds of galax- 
ies (see Wolf et al., in preparation). The redshift quality de- 
pends primarily on apparent magnitude, and is in excellent 
agreement with the simulations of Wolf et al. (2001a). At bright 
limits Mr < 20, the redshifts are accurate to &/(l +z) ~ 0.01, 
and the error is dominated by mismatches between template 
and real galaxy spectra. At the median apparent magnitude 
mR 22, &/(l +z) ~ 0.02. For the faintest galaxies, redshift 
accuracy approaches those achievable using traditional broad- 
band photometric surveys, Sz/(l +z) > 0.05. Restframe colors 
and luminosities are constructed on a galaxy-by-galaxy basis 
by ^-correcting the nearest observed broad band flux. To esti- 
mate this A:-correction, we convolve the best-fit spectrum with 
the filter curves for the Johnson U, B, and V passbands. Typ- 
ical rest-frame color accuracy is — 0.1 mag, corresponding to 
the redshift uncertainty &/( 1 + z) 0.02. Absolute magnitudes 
have an additional error of ~ 0.1 mag (z > 0.5) and ~ 0.2 mag 
(z ^ 0.3) owing to distance uncertainties. 



^This classifier uses a different template set from W03, who use empirical 
templates from Kinney et al. (1996). This completely model-based (somewhat 
unphysical) template set results in superior photometric redshifts to the tem- 
plate set used by W03. Owing to the almost complete age/metalUcity/dust de- 
generacy, we attach no significance to the physical parameters used to derive 
the new template set. 
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Fig . 1 . — The rest-frame U-V color of ~25000 galaxies against the absolute magnitude in V-band, My - 5 log jq /i. We show the distribution of galaxies in nine 
different redshift bins. A fit to the color-magnitude relation of red-sequence galaxies with a fixed slope of -0.08 is shown by the solid line, and the 'Butcher-Oemler' 
style-cut between red and blue galaxies is shown by the dashed line parallel to the early-type galaxy CMR (see §5). The sloping cutoff in the distribution of galaxies 
at the faint end is due to the niR < 24 magnitude limit of the survey. This cutoff cannot be exactly depicted owing to the color dependence in completeness and 
the varying transformation between observed i?-band magnitude and rest-frame V-band magnitude; nevertheless, the dotted grey line shows very schematically the 
rough completeness hmit. The hnes with crosses show the colors and magnitudes of model galaxies with truncated SFHs at constant stellar mass, described in more 
detail in §6.4. Representative error bars are shown. For reference, we show a reddening vector from Calzetti et al. (2000) assuming Ay = 1 mag; a Milky-Way 
extinction curve gives the same vector. Shown in the inset panels is the color distribution of the red peak at each redshift (with the slope of the CMR taken out); 
the dashed line indicates the position of the biweight mean tiiat we adopt as the CMR ridge-line, and the dotted lines show the biweight sigma that we adopt as the 
CMR red sequence's scatter. 
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2.3. Completeness 

We include only galaxies with successful photometric red- 
shift and spectral classifications. Therefore, when construct- 
ing luminosity functions and densities, we must account for 
the sample completeness. This is done using extensive Monte- 
Carlo simulations of the galaxy detection, photometry and clas- 
sification process (WOl); example completeness maps are shown 
in Fig. 7 of W03. For the red galaxies which are the primary 
aim of this paper, the 90% and 50% completeness hmits are at 
roughly ot« 23.2 and 23.8, respectively. 

3 . THE COLOR DISTRIBUTION OF GALAXIES WITH 0.2 < Z < 1 . 1 

One of the main aims of this paper is to explore and discuss 
the color distribution of galaxies in the interval 0.2 < z < 1.1. 
We choose to quantify this in terms of U -V color as a func- 
tion of V-band absolute magnitude, because the U -V color 
straddles the 4000A break and is therefore particularly sensi- 
tive to age and metallicity variations of the stellar populations in 
galaxies, and for consistency with well-known works from the 
literature (e.g., Sandage & Visvanathan 1978; Bower, Lucey, 
& ElUs 1992; Schweizer & Seitzer 1992). In Fig. 1, we show 
the distribution of rest-frame U -V galaxy color against the V- 
band absolute magnitude in nine different redshift bins covering 
0.2<z< 1.1. 

Fig. 1 shows one of the key observational results of this pa- 
per, that the distribution of galaxies in the color-magnitude di- 
agram is bimodal at all redshifts out to z ^ 1 . ^ This confirms 
and extends 9 Gyr back in time the Sloan Digital Sky Survey 
(SDSS) results from e.g., Strateva et al. (2001), Hogg et al. 
(2002), or Blanton et al. (2003a), who established and char- 
acterized the bimodaUty of the color distribution of galaxies in 
the local Universe. Furthermore, this confirms the suggestion 
of Im et al. (2002), who saw hints of a bimodal color distribu- 
tion out to z ~ 1 For immediate low-redshift comparison, we 
show a synthesized U-V color-magnitude diagram for galaxies 
from the SDSS Early Data Release in the Appendix. We defer a 
fuller analysis of this bimodaUty to a future work, but do return 
to this issue briefly in §6.5. 

Beyond the bimodality, it is important to notice two features 
of Fig. 1 . Firstly, the rest-frame U-V colors of blue galaxies 
become redder (at a given absolute magnitude) A(U-V) ^ 0.5 
from z 1 to the present day, probably indicating older mean 
ages and/or larger dust contents, perhaps with a small contri- 
bution from metallicity evolution at a given luminosity. In ad- 
dition, the number density of blue galaxies, especially at high 
luminosities, significantly evolves. Because of the different vol- 
umes sampled by COMBO- 17 in each redshift interval, this 
strong evolution in the number density of luminous blue galax- 
ies is not apparent in Fig. 1. W03 studied this issue in more 
detail; when volume-corrected, the number density of faint star- 
forming galaxies (spectral Type 4) remains almost unchanged, 
the abundance of luminous galaxies with starburst spectra drops 
precipitously from redshift uiuty to the present epoch (Fig. 16 
of W03; also, e.g., Cowie et al. 1996). This result, which is dis- 
cussed much more by W03, also agrees with the CNOC2 (Lin 

'it is worth noting that the full range of rest-frame U-V colors is covered 
smoothly by the template spectra, making it unlikely that the algorithm arti- 
ficially focuses galaxies away from the gap between the red and blue galaxy 
sequences. Furthermore, red-sequence galaxies at z ~ 1 have (/- V ~ 1.2, and 
are the same color as gap galaxies at z ~ 0.2, arguing further against artificial 
de-focusing of tZ-V ~ 1.2 galaxies to redder or bluer colors. 

"•The Deep2 survey with Keck also appears to find this bimodality (B. 
Weiner, 2003, priv. comm.) 
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Fig. 2. — Color evolution of the 'red sequence', represented by the intercept 
of the CMR fits at My - 5 log | (, /; = -20 (solid circles) as a function of redshift. 
The error bars show the biweight scatter in each redshift bin. These highly 
conservative errors are much larger than the formal errors, and likely overesti- 
mate the systematic errors from calibration, template/galaxy SED mismatches 
and small redshift focusing issues, which we expect to be less than 0.1 mag in 
U-V color. The lines show the expected color evolution of single-age stel- 
lar populations with different formation redshifts and metallicities, as given in 
the Figure legend. The diamond with error bars (offset from zero redshift for 
clarity) shows Aie z = U-V CMR zero point synthesized from SDSS ugr 
data analyzed in the same way as these data (see the Appendix for details). 
The naked error bars show the U-V CMR zero point for the Abell 754 galaxy 
cluster (Mcintosh, Rix, & Caldwell 2004), which is also consistent with the 
Nearby Field Galaxy Survey (Jansen et al. 2(X)0) early-type galaxy CMR. The 
black naked error bar shows the AbeU 754 scatter around the CMR, and the 
expected systematic uncertainty is shown by a grey error bar. The grey line 
shows a fit to aU of the data (including the Abell 754 and SDSS colors). 



et al. 1999) and CFRS (Lilly et al. 1995) surveys, who found 
sttong evolution in the abundance of star-forming, blue galax- 
ies with redshift. We choose not to elaborate further on the 
evolving colors or luminosities of star-forming galaxies at this 
stage. 

Secondly, bearing in mind the > 0. 1 mag errors in ^-corrected 
rest-frame U-V colors, the red galaxies form a relatively well- 
defined sequence in the color-magnitude plane for 0.2 < z < 
1.1. Furthermore, the mean color of this red sequence evolves 
with redshift, in the sense that galaxies become bluer at a fixed 
luminosity. This evolution will be quantified in §4. Both ob- 
servations agree with studies of galaxies in clusters at redshifts 
up to unity (e.g.. Bower, Lucey, & Ellis 1992; Kodama & Ari- 
moto 1997; Terlevich, CaldweU, & Bower 2001; van Dokkum 
& Franx 2001), but we extend these results by showing that 
they apply to volumes which average over many different envi- 
ronments. The red sequence has a significant slope in the color- 
magnitude plane, therefore can be rightfully termed a color- 
magnitude relation (CMR) in all redshift bins up to at least 
z ~ 0.8. At z > 0.8, the magnitude range is so narrow that a 
constant color is also quite plausible. It is the CMR of the red 
sequence galaxies that we study in detail in the rest of this pa- 
per. 
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4. THE COLOR-MAGNITUDE RELATION OF RED GALAXIES 

We now quantify the redshift evolution of the CMR by exam- 
ining its zero point. In order to ensure consistency between the 
zero point estimates for different bins, we constrain the slope 
of the CMR to the locally-determined value for cluster galax- 
ies d(U-V)/dMv = -0.08 (Bower, Lucey, & Elhs 1992; Ter- 
levich, Caldwell, & Bower 2001). This slope appears com- 
patible with the CMRs presented in Fig. 1 (a quantitative as- 
sessment of this statement is challenging without developing a 
robust description of the entire bimodal distribution, which is 
beyond the scope of this work; see, e.g., Baldry et al. 2004), 
and indeed is compatible with the slope of -0.05 derived by 
Schweizer & Seitzer (1992) for E and SO galaxies in field envi- 
ronments. We determine the robust biweight mean color, with 
the locally-determined CMR slope subtracted, of all galaxies 
with U-V > 1.0 (independent of redshift, so as to roughly iso- 
late the red sequence). We illustrate the procedure by show- 
ing the CMR-subtracted color histogram of galaxies in the inset 
panels of Fig. 1 . The biweight mean and sigma of the unbinned 
sample are denoted by the dashed and dotted lines respectively. 
It becomes clear that the peak of the distribution can be de- 
termined to within a few hundredths of a magnitude, which is 
sufficient for our purposes. 

The value of our fit to the CMR color at Mv - 5 logj,, h = -20 
is shown as a function of redshift in Fig. 2 and is tabulated in 
Table 1. It is clear that this population evolves considerably 
between z = I.I and 0.2; the galaxies become redder by ~ 0.3 
mag. For completeness, a linear fit to the evolution (including 
the local comparison points) is {U -V)mv=-2() = 1.40-0.31z, 
as shown by the solid grey line. Some idea of the systematic 
uncertainties is given by a fit to the COMBO- 1 7 points only: 
([/-y)Mv=-20 = I.48-0.40z (the fit is not shown for clarity). 
The biweight RMS for the red sequence is given by the error 
bar. Overplotted are the evolving U-V colors expected for 
single-age populations with differing formation redshifts. We 
use the Pegase stellar population synthesis model with subso- 
lar metalhcities to predict colors at a constant absolute magni- 
tude. Galaxies with My - 5 logjQ h = -20 at redshift zero would 
be brighter at redshift one, owing to passive evolution of the 
stellar population. Therefore, we apply a small color correc- 
tion to the passive evolution models to account for the fact that 
the galaxies with My -Slog^^h = -20 at redshift one are from 
a different, bluer part of the CMR. The evolution of the model 
colors are insensitive to the choice of stellar initial mass func- 
tion (IMF), for which we adopt the Salpeter (1955) parameter- 
ization. We show also the color of the CMR intercept from 
two local samples. The color intercept of morphologically- 
classified early-type galaxies in the Abell 754 galaxy cluster 
is shown by the naked error bars (Mcintosh, Rix, & Caldwell 
2004). This color intercept is consistent with early-type galax- 
ies from a wide range of environments from the Nearby Field 
Galaxy Survey (Jansen et al. 2000; Mcintosh, Rix, & Caldwell 
2004). The U-V CMR intercept synthesized from SDSS Early 
Data Release (EDR; Stoughton et al. 2002) u, g, and r-band data 
is shown by the diamond with error bars, which is offset from 
zero redshift for clarity. The SDSS EDR sample is discussed in 
the Appendix. 

The color evolution of the COMBO- 17 CMR intercept (the 
(U-V) at My - 5 log JO ^ = -20) is consistent with the expecta- 
tions of passive evolution of ancient stellar populations. How- 
ever, the COMBO- 17 points do seem slightly high compared 
to the local data at z < 0.5. Absolute caUbration of surveys is 



challenging at even the ^ 10% level, especially as all of the 
data shown in Fig. 2 have been transformed to U-V from other 
passbands, except for the Abell 754 data from Mcintosh, Rix, 
& Caldwell (2004). Furthermore, the number of red-sequence 
galaxies is low at z < 0.4, increasing the uncertainties consider- 
ably (see Fig. 1). Finally, and perhaps most importantly, color 
gradients will affect this analysis at some level. The COMBO- 
17 colors are measured using adaptive-aperture photometry, 
and are essentially small-aperture central colors. For the most 
distant galaxies this will approximate total color For the near- 
est COMBO- 17 galaxies at z < 0.4, we may sample just their 
inner, redder parts. Taking typical color gradients from Peletier 
et al. (1990) as a guide, this could easily lead to a ^ 0.1 mag 
reddening of the CMR intercept at low redshift compared to 
the passive evolution expectation from high redshift. The lo- 
cal surveys we show use colors derived for between 50% and 
100% of the galaxy light, and this would also lead to an offset 
between the lowest-z COMBO- 17 and local color intercepts at 
around the ~ 0.1 mag level. Both of these effects are arguably 
seen in the data. Nevertheless, taken as an ensemble, the CMR 
intercepts do appear to redden with time in a way that is con- 
sistent with the expectations of passive evolution. Yet, even 
extending large samples to z ^ 1, it is not possible to constrain 
meaningfully the formation redshift of the stars in early-type 
galaxies using this type of information, owing both to the age- 
metallicity degeneracy, and to systematic model uncertainties 
of at least ~ 10% in the colors of old stellar populations (Char- 
lot, Worthey, & Bressan 1996). This is illustrated in Fig. 2, 
where we choose stellar populations of different ages, with their 
metallicity at My - 5 log jq h = -20 adjusted to keep the average 
U-V color evolution roughly constant. 

It is worth briefly discussing the scatter in the CMR of early- 
type galaxies. In the local Universe, the scatter in the CMR 
is observed to be small (< 0.04 mag; Bower, Lucey, & Elhs 
1992; Terlevich, Caldwell, & Bower 2001, alfliough Schweizer 
& Seitzer 1992 find a scatter of ^0.1 mag in a more field 
galaxy-dominated sample). This has been interpreted as evi- 
dence for a small age spread in local Universe early-type galax- 
ies, although Trager et al. (2000) argue that age and metalhc- 
ity effects can conspire to produce a small CMR scatter even 
with large, anti-correlated scatter in age and metallicity. We 
find rather larger scatter (Fig. 1), mostly because our photomet- 
ric redshift errors of Az ~ 0.03 translate into errors in the de- 
rived rest- frame colors of A(U-V) ~ 0.1 mag. In addition, the 
red sequence at z < 0.4 is poorly defined, leading to larger un- 
certainties. Owing to these, and other systematic uncertainties 
(such as calibration and mismatches between real galaxy spec- 
tra and the template spectra used to derive the A:-correction), we 
show the full zLla scatter in Fig. 2. In any case, these sources 
of uncertainty do not significantly affect our analysis because 
the color selection cut of A{U -V) = 0.25 mag that we choose 
later in §5 is larger than the scatter at all redshifts (< 0.2 mag). 

5. THE LUMINOSITY DENSITY OF RED-SEQUENCE GALAXIES 

Because the stellar populations of red-sequence/early-type 
galaxies may have formed much earlier than they were assem- 
bled into recognizable early-types (Baugh, Cole, & Frenk 1996; 
van Dokkum & Franx 2001), the evolution of the CMR, or of 
the stellar mass-to-light ratios of red galaxies, is not a sensitive 
diagnostic of when these galaxies attained their present con- 
figuration. A much more sensitive diagnostic is the density or 
luminosity functions of red-sequence galaxies at redshifts less 
than unity. 
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Fig . 3 . — Redshift evolution in tlie rest-frame B-band luminosity function of red-sequence galaxies. Black solid lines with data points show the Schechter Function 
fits and y/Vmax luminosity functions of COMBO- 17 red-sequence galaxies in each redshift interval. The dotted line shows the Schechter fit to the low-redshift SDSS 
comparison sample (see the Appendix for more details). A single faint-end slope a = —0.6 is derived by a single STY fit to the COMBO-17 data at all redshifts. 
When deriving the if>* estimate, the three least luminous bins are ignored. The arrow shows the magnitude cutoff used for the analysis of the luminous red-sequence 
galaxy population in §6.4, assuming zf = 3. 



In §§3 and 4, we have seen that it is possible to define the 
early-type galaxy population in an empirical, model-independ- 
ent fashion by exploiting the bimodality of the galaxy color dis- 
tribution out to z ~ 1 (Fig. 1). We therefore choose to study 
galaxies around the red galaxy CMR, defining these as the red- 
sequence galaxy population at any given redshift. More quan- 
titatively, we define red-sequence galaxies as being redder than 
(U-V) = 1.15-0.3U-0.08(Mv-51ogio/i-i-20) (the dashed 
Une in Fig. 1). This line is 0.25 mag bluewards of the linear 
fit to the evolution of the mean color of a My - 5 logjg h = -20 
galaxy from Fig. 2 (the solid grey line)^. We choose this def- 
inition merely to have a criterion that is a continuous function 
of redshift. Using a discontinuous color-cut which is 0.25 mag 
bluewards of the empirical CMR zero-point defined within each 

^This fit includes the data points also. A color cut derived using just 
the COMBO-17 data (which therefore includes the small bias caused by color 
gradients) is: {U-V) = 1.23 -0.40z -0.08(Mv -5 logj,,/? 4-20). Adoption of 
this color cut in what follows does not change the results to within the errors. 



Az = 0.1 bin yields results thai agree within the errors. Our 
color-cut philosophy is similar to that of Butcher & Oemler 
(1984), except that we choose to keep the red galaxies rather 
than the blue ones. The advantages of this definition of early- 
type galaxies are many. Firstiy, it is straightforward to match 
this selection criterion in theoretical models. For example, even 
if the model CMR is curved, one can simply go ^0.25 mag red- 
wards of the curved locus in U-V space and choose galaxies 
redder than that curved locus (although as long as the 'gap' 
between the red and blue peaks is sparsely populated, like in 
the observational data, the exact choice of color cut makes little 
difference as long as the cut lies in the 'gap'). Secondly, it auto- 
matically accounts for evolution in a model-free way by fitting 
the CMR ridge-Une. Thirdly, this definition has a well-defined 
physical meaning (galaxies that have colors consistent with old 
stellar populations at a given redshift). A final advantage is 
this definition's insensitivity to differences in photometric zero 
points. For example, the local SDSS comparison sample and 
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Table 1 

Red-Sequence Galaxy Rest-frame B-band Luminosity Function Fits and CMR Intercept 



(z) 


logh 


(/}* X 10""* 


IB X lO'LfTi 


J obs 








(Vega mag) 


(/i/Mpc)-^ 


(h/Mpc^) 








0.25 


-19.94±0.22 


71±16, 20 


8.54±1.88, 2.34 


1.0 


-0.297 


1.41 ±0.22 


0.35 


-19.98±0.19 


56±23, 12 


6.87±2.81, 1.42 


1.0 


-0.291 


1.35±0.19 


0.45 


-19.92±0.17 


62±9, 10 


7.27±1.10, 1.23 


0.99 


-0.300 


1.30±0.18 


0.55 


-20.25±0.17 


44±2, 6 


6.94±0.38, 1.00 


0.99 


-0.251 


1.24±0.16 


0.65 


-20.34±0.17 


35±7, 5 


6.09±1.14, 0.78 


0.97 


-0.315 


1.16±0.12 


0.75 


-20.44±0.16 


43±3, 5 


8.16±0.54, 0.94 


0.94 


-0.456 


1.14±0.14 


0.85 


-20.33±0.14 


50±11,5 


8.43±1.91,0.89 


0.88 


-0.675 


1.15±0.15 


0.95 


-20.64±0.17 


20±5, 2 


4.50±1.11,0.44 


0.85 


-0.767 


1.10±0.13 


1.05 


-20.74±0.22 


16±3, 1 


3.87±0.64, 0.35 


0.77 


-1.019 


1.09±0.12 



Note. — STY fits in 9 redshift intervals where the faint end slope a = -0.6 ±0.1 was fit to the entire sample. Listed are M*, (j)* with 
two estimates of the cosmic variance uncertainty; from field-to-field variation divided by -\/3, and following the prescriptions of 
Somerville et al. (2004). The luminosity density Jb is given, again with cosmic variance uncertainties, and further hsted is the 
fraction of the luminosity density that is constrained by the measurements, and the covariance between </>* and M* . Also given is the 
intercept of the CMR at {U -V) color at My = -20 and the robust biweight RMS scatter. 



COMBO- 17 sample may have offset U -V color scales (Fig. 
2); yet, the B-band luminosity density of galaxies on the CMR 
is insensitive to this small color offset, meaning that the lumi- 
nosity density of the local SDSS sample can be meaningfully 
compared with the rest-frame B-band luminosity density of the 
COMBO- 17 dataset. We note that this definition contrasts with 
W03, who adopted spectral type definitions based on present- 
day spectral energy distribution types. With this new approach, 
we are taking into account in a natural and empirical way the 
effects of evolution in red sequence galaxy color. 

Defining galaxy types in this way, we now explore the evolu- 
tion of the luminosity function of 4690 0.2 < z < 1 . 1 COMBO- 
17 red-sequence galaxies in Fig. 3. We estimate the galaxy 
luminosity function in different redshift bins following W03 
using the V/Vmax formalism (e.g., Felten 1977), using our de- 
tailed estimates of galaxy completeness as a function of sur- 
face brightness, magnitude, redshift and spectral type (see, e.g.. 
Figs. 6 & 7 of W03). Error bars on the V/Vmax data points 
take into account Poisson errors only (i.e., not accounting for 
redshift, spectral type and magnitude errors). Schechter (1976) 
functions were fit to the red-sequence galaxy sample using the 
method of Sandage, Tammann, & Yahil (1979, STY hereafter).^ 
We choose to fit the entire red-sequence galaxy sample at all 
redshifts with one common value of the faint-end slope a; in 
practice, a is determined mainly by the lower redshift galaxy 
samples and is -0.6 ±0.1 in this case, in excellent agreement 
with local observations (e.g., Madgwick et al. 2002; Bell et al. 
2003). The magnitude at the 'knee' of the luminosity function 
(M*) is redshift, but not field- to-field, dependent, and the den- 
sity normaUzation (^* is fit separately to each field at each red- 

*In studies where the luminosity function is derived from a redshift survey 
with considerable 'thickness', such that galaxies populating the faint end come 
from different volumes as the galaxies populating the bright end, the V/Vmax 
can give misleading results owing to large scale structure. The STY estimator 
does not suffer from this limitation (Sandage, Tammann, & Yahil 1979). In 
this paper, we split the COMBO- 17 survey data into many thin shells, therefore 
the faint and bright galaxies in each redshift interval come from nearly iden- 
tical volumes. Thus, the V/Vmax luminosity function estimates should also be 
relatively robust to large scale structure. 



shift (see Table 1 for fit parameters). The Schechter luminosity 
function to the local comparison sample is shown by a dotted 
line (see the Appendix for more details). 

W03 found strong evolution of 'early-type galaxies', de- 
fined via redshift-independent spectral energy distribution cri- 
teria (which did not take into account galaxy evolution). In 
contrast, one can see from Fig. 3 that the luminosity function 
evolution of early-type galaxies defined in this way is relatively 
mild. The characteristic magnitude M* seems to brighten and 
the number density (jf seems to decrease somewhat with in- 
creasing redshift. 

We quantify this evolution in Fig. 4, where we show the red- 
shift evolution of the red-sequence galaxy rest-frame B-band 
luminosity function parameters. In the left and center pan- 
els, we show the evolution of the magnitude of the 'knee' of 
the luminosity function, M%, and the density normalization of 
the luminosity function, cf)* . The solid circle denotes the low- 
redshift SDSS EDR comparison sample, selected using synthe- 
sized U — V colors (Bell et al. 2003, see the Appendix), and k- 
corrected and evolution-corrected to redshift zero. The smooth 
grey fines show the expectation of passive evolution models, 
normalized to pass through the COMBO-17 Mg evolution. The 
model evolution is insensitive to stellar IMF at the 0.1 mag 
level. Error bars in show the 68% confidence interval. The 
error bars for (j)* attempt to show the expected effect of cosmic 
variance, and were derived using two methods. The black error 
bars show the field-to-field variations divided by -\/3, whereas 
the grey error bars show predictions from the prescription of 
Somerville et al. (2004), which over this redshift range for our 
sample of galaxies depends primarily on the comoving volume. 
As Fig. 4 and Table 1 show, these totally independent error es- 
timates agree quite well, demonstrating that we have accounted 
for the effects of cosmic variance on our results. It is worth 
noting that M% and </>* estimates are tightly coupled to give an 
invariant result for the total luminosity density Jb ^. One can 

^This is a possible explanation for the offsets of the SDSS z, = data point 
from extrapolations from the COMBO-17 data, as the SDSS control sample 
may have a sUghtly different shape, giving low Mg and high 0* estimates. 




The Evolution of Early-Type Galaxies 
-1.8 

-2.0^ 

-2,2 f i^To 



\ -2.4 
«- 

^ -2.6 
o 

-2.8 

-3.0 



C) 



C) 



C) 



C). 




0.0 0.2 0.4 0.6 0.8 1.0 1 
z 



.2 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
z 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
z 



Fig. 4. — The evolution of the rest-frame B-band luminosity function of the red-sequence galaxy population from COMBO-17. Left: The evolution of Mg, the 
'knee' of the luminosity function. We show the expectations of a passive evolution model as smooth grey lines (with formation redshift shown in the Figure legend). 
The soUd point shows the Mg of the SDSS low redshift comparison sample, ttansformed from ugr data (see the Appendix for more details). Center: The evolution 
of (p*: again, open circles denote COMBO-17 data and the solid point the local SDSS sample. The error bars attempt to indicate uncertainty owing to cosmic 
variance; the black error bars show the observed field-to-field variation divided by \/3, whereas the grey error bars show the predicted cosmic variance using the 
prescription of SomerviUe et al. (2004). A passive evolution model predicts no evolution in the number density with redshift (grey line). Right: The rest-frame 
B-band luminosity density per comoving Mpc^ is shown by open circles, with error bars again denoting the uncertainties from cosmic variance. To reduce the 
cosmic variance uncertainty, we show also the averaged 0.2 < z < 0.9 js by the soUd small black circle, which has only a 13% uncertainty from cosmic variance. 
The upwards-pointing triangles denote a lower limit to the luminosity density, where only the observed luminosity density is accounted for (i.e., no extrapolation to 
zero luminosity). The filled grey point denotes the luminosity density of spectral-type selected early-type galaxies in the 2dFGRS. We show also the blind prediction 
of the Cole et al. (2000) semi-analytic galaxy formation model as a thick soUd grey curve, where model galaxies were selected using exactiy the same color-cut 
methodology as used in this paper 



immediately see that the Mg evolves, like the colors, in a way 
that is consistent with the evolution of an ancient stellar popula- 
tion. In contrast, the 0* seems to strongly evolve in the interval 
< z < 1, in clear disagreement with passively-evolving or pure 
luminosity evolution models. 

Much more robustly measured is the total luminosity density 
of the summed red-sequence galaxy population. The luminosity 
density is estimated by integrating the best-fit Schechter Func- 
tion to infinity using j = (j)*L*T{2 + a). We show the average 
luminosity density at each redshift from red-sequence galaxies 
(open circles), again including an error bar showing the uncer- 
tainty from cosmic variance. We show also a lower limit to the 
luminosity density, derived only from the observed luminos- 
ity density without extrapolating to zero luminosity (upwards- 
pointing triangles). We show also the SDSS (solid circle) lu- 
minosity density, and the 2dFGRS summed luminosity density 
of spectral-type selected early-type galaxies (Madgwick et al. 
2002) as a filled grey circle. The smooth grey lines show the 
expectation of passive evolution models as shown above. 

This figure shows a few key results. Firstly, it is clear that 
cosmic variance is a considerable source of uncertainty, even 
for a survey with 3 independent fields of ~ 0.25 square de- 
gree each. Cosmic variance is the Umiting uncertainty for our 
and similar studies of the < z < 1 luminosity function evo- 
lution of red-sequence or early-type galaxies. Furthermore, 
SomerviUe et al. (2004) shows that more than an order of mag- 
nitude increase in comoving volume is required to reduce the 
cosmic variance uncertainties by a factor of two, if the areas 
imaged are contiguous. Alternatively, a factor of fotir more 
widely-separated areas could give the same decrease in un- 
certainty. Therefore, it will be necessary to cover substan- 
tially larger areas to similar depths to sigiuficantiy improve 
on these results. Secondly, and most importantly, the rest- 
frame B-band luminosity density in red-sequence galaxies re- 
mains more or less constant over the interval < z < 1. In 
an effort to reduce the effects of cosmic variance, we estimate 
the average luminosity density in red-sequence galaxies in the 



interval 0.2 < z < 0.9, at the risk of over-binning the data: 
jB(0.2 < z < 0.9) = 7.6 ± 1.0 X lO^Lo/zMpc^l The observed 
luminosity density falls short of the passive evolution predic- 
tion by at least a factor of two by z ~ 1, although cosmic vari- 
ance and uncertainty in the passive evolution models preclude 
a more accurate assessment. 

Of some importance in establishing this result is the local 
estimate of red sequence luminosity density. This was the pri- 
mary motivation for carefully constructing a SDSS comparison 
sample. As discussed in Bell et al. (2003), this estimate is accu- 
rate to 10% in a systematic and random sense, owing largely to 
the excellent completeness properties of SDSS (used for galaxy 
selection) and the Two Micron AH Sky Survey (which is used 
to normahze the luminosity function to the all-sky number den- 
sity of \0 < K < 13.5 galaxies). Furthermore, the U -V col- 
ors were synthesized to allow identical selection criteria to be 
appUed to the low-redshift and COMBO-17 samples, strongly 
limiting that source of systematic uncertainty. Furthermore, the 
2dFGRS estimate of Madgwick et al. (2002) agrees with our 
SDSS estimate, lending further credibility to the low-redshift 
comparison data. Therefore, we believe that it is safe to con- 
clude that the rest-frame fi-band luminosity density of color- 
selected early-type galaxies does not significantly evolve in the 
interval < z < 1 . Furthermore, because a passively-evolving 
population would have to fade by a factor of two or three from 
redshift one to the present day in rest-frame B-band, this non- 
evolving Jb corresponds to an increase in stellar mass in red 
galaxies of a factor of two to three from redshift unity to the 
present day. 

6. DISCUSSION 

6.1. Comparison with previous work 

The data presented here are based on a larger galaxy sam- 
ple than all previous work in this field combined. Further- 
more, we adopt an empirically-motivated definition of 'early- 
type galaxy' which accounts for evolution of the stars in these 
galaxies in a natural and model-independent fashion. Never- 
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theless, on the whole there is pleasing agreement between our 
conclusions and the results of many of the previous studies 
of this topic. The essentially passive evolution of the stars 
in early-type galaxies (not the galaxy population as a whole), 
as estimated from the color intercept of the CMR, is in excel- 
lent agreement with the colors (e.g., Kodama & Arimoto 1997; 
Kodama, Bower, & Bell 1999; van Dokkum et al. 2001), line 
strengths (e.g., Kuntschner 2000; Kelson et al. 2001; Ziegler et 
al. 2001), and stellar M/Ls (e.g., van Dokkum & Franx 2001; 
Treu et al. 2002) of individual early-type galaxies between 
< z < 1, as extensively discussed in the introduction. 

There is also good agreement between our detection of a fac- 
tor of two evolution in stellar mass on the red sequence and 
other works in the literature. Morphologically-selected surveys 
(Im et al. 1996; Schade et al. 1999; Menanteau et al. 1999; 
Brinchmann & Elhs 2000; Im et al. 2002) found little evidence 
for density evolution given the small sample sizes (< 150 galax- 
ies) and cosmic variance, and could not rule out passive or mild 
evolution of the type that we see. Noting that the most luminous 
galaxies are red at all redshifts up to unity (Fig. 1), our results 
also agree with a factor of a few decrease in the number of very 
luminous galaxies in A'-band selected samples out to z ~ 1 (e.g., 
Drory et al. 2001). 

Most color-selected surveys also agree with our conclusions. 
Lilly et al. (1995) study evolution in the luminosity function 
of galaxies with colors redder than a present-day Sbc galaxy. 
They find very little evolution, with large uncertainties from 
small number statistics and large-scale structure, which is con- 
sistent with both our result and passive evolution. Kauffmann et 
al. (1996) re-analyse this result, finding a factor of three evolu- 
tion in galaxies with the colors of early-types, although Totani 
& Yoshii (1998) point out that spectroscopic incompleteness 
affects this result. Lin et al. (1999) agree with this picture us- 
ing a similar technique with a larger sample of galaxies over 
a more restricted redshift range. Pozzetti et al. (2003) used 
the K2Q survey to estimate the evolution of red galaxies out 
to z ~ 1 ; using a total galaxy sample of only 546 galaxies from 
52 square arcmins, they find mild evolution in the number of 
luminous red galaxies, but with error bars of more than a fac- 
tor of two (not accounting for the dominant uncertainties from 
large scale structure). Cimatti et al. (2003) study the properties 
of extremely red objects (EROs) in K2Q (galaxies with the col- 
ors of passively evolving populations), finding 2-Ax less EROs 
than predicted from passive models (although with large cos- 
mic variance uncertainty); furthermore, many of these galaxies 
have disturbed or disk-like morphologies. Chen et al. (2003) 
have used the Las Campanas Infrared Survey to estimate the 
luminosity density evolution of red galaxies in the rest frame 
/?-band, finding at most a factor of a few evolution in red galaxy 
density, although they were not able to rule out passive evolu- 
tion with their sample. We are consistent with these results, but 
have been able to place much stronger constraints based on our 
factor of at least three improvement in sample size and factor 
of two larger sky coverage, compared to the largest of the pre- 
vious works. Importantly, these previous works often neglected 
the uncertainties stemming from large-scale structure; whereas, 
we have quantified and accounted for this dominant source of 
uncertainty. 

Finally, we should compare with our own earlier work, W03, 
where we explored the same COMBO- 17 dataset, but used a 
concepually different, non-evolving definition of 'early-type 
galaxies' (their type 1). Their non-evolving spectral typing 



corresponds to constant rest-frame color cuts at all redshifts; 
their type 1 corresponds to fZ-V > 1 .35, and their type 2 corre- 
sponds to 1.35 > U — V > 0.95. Thus, type 1 galaxies have the 
colors of present-day luminous elliptical and lenticular galax- 
ies, but are significantly redder than the reddest galaxies at 
higher redshift. Consequently, W03 found rapid evolution of 
type Is, in the sense that galaxies with the colors of present-day 
early-types are extremely rare in the distant Universe (but pri- 
marily because of passive evolution in their colors). Because of 
these definition differences, the statement by W03 that galaxies 
redder than a fixed threshold have a rapidly evolving luminos- 
ity density is consistent with our statement of a roughly non- 
evolving rest-frame B-band luminosity density in red-sequence 
galaxies. 

In summary, all of our results are consistent with those pre- 
sented to date in the hterature, but with much greater precision 
and understanding of the uncertainties owing to our factor of 
> 3 larger sample and sky coverage. 

6.2. The Nature of Red Galaxies 

We now turn to the nature of galaxies in this red population: 
are these galaxies red because they contain almost exclusively 
old stars, or is there a significant contribution from more dusty, 
star-forming galaxies? 

At low redshift, the red population is overwhelmingly com- 
posed of morphologically early-type galaxies (e.g., Sandage & 
Visvanathan 1978; Bower, Lucey, & Elhs 1992; Schweizer & 
Seitzer 1992; Terlevich, Caldwell, & Bower 2001; Hogg et al. 
2002). For the SDSS survey, Strateva et al. (2001) find that 
~80% of the red galaxies are earlier than Sa in morphological 
type. 

At intermediate redshift, most studies have targeted galaxy 
clusters where most red or non-star-forming galaxies are spheroid 
dominated (e.g.. Couch et al. 1998; van Dokkum et al. 2000; 
van Dokkum & Franx 2001). Furthermore,most of the lumi- 
nosity density of morphologically-classified early-type galaxies 
in the Hubble Deep Field North at z ~ 0.9 is in red-sequence 
galaxies (Kodama, Bower, & Bell 1999). Recently, Bell et 
al. (2004) used the GEMS (Galaxy Evolution from Morphol- 
ogy and SEDs Rix et al. 2004) HST Advanced Camera for 
Surveys dataset to explore the morphology of roughly 500 
0.65 > z > 0.75 red-sequence galaxies, finding that > 80% 
of the luminosity density in z ~ 0.7 red-sequence galaxies is 
from morphologically-classified early-type (E, SO and some Sa) 
galaxies. 

The nature of red galaxies at higher redshifts 1 ^ z < 2 is 
somewhat less clear. A number of studies have explored the 
nature of 'extremely red objects' (EROs), which have very red 
colors consistent with passively-evolving ancient stellar popu- 
lations with 1 < Z ^ 2 (but can in principle be a mixed bag from 
a variety of different redshifts). Cimatti et al. (2002) use deep 
spectroscopy to split EROs roughly 50:50 between passively- 
evolving red-sequence galaxies and dusty star-forming galaxies 
in a A" < 19.2 sample of 45 EROs. Complementary analyses are 
presented by Yan & Thompson (2003), Moustakas et al. (2004), 
Gilbank et al. (2003), and Cimatti et al. (2003), who demon- 
strate that there is roughly an equal three-way split between the 
Hubble Space Telescope (HST) morphologies of EROs, such 
that ~ 1/3 are bulge-dominated, ^ 1/3 are disk-dominated, and 
^1/3 have disturbed morphologies. Thus, luminosity densities 
derived from red color-selected 1 ^ z < 2 samples may overes- 
timate the stellar mass in early-type galaxies. 
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This evidence suggests that most red galaxies to z < 1 are 
morphologically early-type with dominant old stellar popula- 
tions. In contrast, a significant fraction of distant z > 1 red 
galaxies may have red colors not because of a dominant old 
stellar population, but rather because of dust (i.e., edge-on disks 
or dusty starbursts). Bearing this in mind, we caution that our 
results should be interpreted only as upper limits on the abun- 
dance of galaxies dominated by ancient stellar populations until 
the advent of more definitive morphological data and/or spectra 
for sufficient numbers of galaxies with 0.2 < z < 1. 

6.3. The Predicted Luminosity Evolution of Hierarchical 
Model Galaxies 

In §5 we showed that the rest-frame fi-band luminosity den- 
sity of red-sequence galaxies does not significantly evolve since 
z ^ 1. This is at variance with the expectations of a com- 
pletely passive evolution model (where red-sequence galaxies 
fully form at z 1 and simply age to the present day), which 
would predict a factor of two or three decrease in the rest-frame 
B-band luminosity density in red galaxies from redshift one to 
zero. It is interesting to compare our results with the expecta- 
tions of hierarchical models of galaxy formation and evolution, 
to check if these models would predict a roughly non-evolving 
B-band luminosity density since z ~ 1 in red galaxies. 

There are two main predictions of hierarchical models of 
galaxy formation and evolution. Firstly, early-type galaxies in 
the field should have stellar population of a somewhat younger 
mean stellar age, and they should be assembled later than their 
clustered counterparts, because the mergers which give rise to 
the early-types happened at earlier times in clusters than in the 
field (where they continue even to the present day; e.g., Baugh, 
Cole, & Frenk 1996). We have chosen not to explore this issue 
at present, as constructing local density estimators with a red- 
shifl precision of Az ^ 0.03 is non-trivial. Other studies may 
have detected an offset, in the sense that field early-types may 
be younger than cluster early-types, using Une strength (Guz- 
man et al. 1992; Bemardi et al. 1998) or fundamental plane con- 
straints (Treu et al. 2002), although Hogg et al. (2004) see no 
environmental dependence in the CMR in the SDSS at the 0.01 
mag level. Nevertheless, progenitor bias and the age-metalhcity 
degeneracy are significant observational worries in studies of 
this kind, and more work is required to firm up these interest- 
ing results. A particularly exciting avenue of research is study 
of mass-selected galaxies from gravitational lensing studies. In 
these studies, the sample is selected to have high stellar-hdark 
masses, without regard for their photometric properties (or in- 
deed, for their environments). First results from these surveys 
have been controversial but promising. Kochanek et al. (2000) 
and Rusin et al. (2003) find no evidence for younger stellar 
populations in (the same) mass-selected sample, compared to 
galaxy clusters, van de Ven, van Dokkum, & Franx (2003) 
re-analyze these data using a more model-dependent but sensi- 
tive approach, accounting for progenitor bias, finding younger 
luminosity-weighted ages than clustered massive galaxies, in 
agreement with hierarchical model expectations. 

The second main prediction of hierarchical models of galaxy 
formation is a gradual increase in the stellar mass of the early- 
type/red- sequence galaxy population until the present day, be- 
cause of constant build-up in stellar mass from galaxy mergers 
(Baugh, Cole, & Frenk 1996; Kauffmann et al. 1996). In Fig. 
4, we show a blind prediction of evolution of the rest-frame B- 
band luminosity density of red-sequence galaxies from the Cole 



et al. (2000) fiducial model using the same color-magnitude re- 
lation methodology adopted in this paper. The models are con- 
strained to match the local luminosity function plus the mor- 
phological mix of L > L* galaxies in the local Universe, there- 
fore the approximate match between the models and the ob- 
servations at redshift zero is largely by design. Observations 
of luminosity functions at non-zero redshifts were not used at 
any stage to constrain model parameters, therefore the excellent 
agreement between the model and the data out to z ~ 1 is a sig- 
nificant achievement. It is unclear how well other models (e.g., 
Kauffmann et al. 1999; Somerville & Primack 1999) will repro- 
duce these data, owing to the well-known inability of models 
of this type to match both the dynamics and luminosity func- 
tions of galaxies simultaneously. The Cole et al. (2000) model 
tunes to match the luminosity function of galaxies well at the 
present day, and is therefore the suitable model with which to 
compare luminosity density evolution. The Somerville & Pri- 
mack (1999) and Kauffmann et al. (1999) models tune to repro- 
duce galaxy dynamics at the expense of the galaxy luminosity 
function, and are less suitable for comparisons of this type (but 
are more suitable for exploring the evolution of galaxy dynam- 
ics). A more thorough discussion of the stellar mass, dynamical 
and morphological evolution of the red-sequence galaxy popu- 
lation is clearly warranted, but is beyond the scope of this pa- 
per. Nevertheless, the current zero-order comparison is rather 
encouraging; models that are tuned to reproduce the present- 
day distribution of galaxy luminosities predict no evolution in 
the luminosity density of red-sequence galaxies, in agreement 
with our observations. 

6.4. The Origin of Luminous Red-Sequence Galaxies 

What might these data tell us about the origin of luminous 
red-sequence galaxies? An intriguing observation is that at all 
redshifts, there are very few blue galaxies luminous enough to 
fade into the brightest red-sequence galaxies (Fig. 1). In this 
section, we explore what this may imply. 

As a useful tool for this discussion, we construct a sequence 
of previously star-forming Pegase model galaxies who have 
had their star formation suddenly stopped at some point in the 
past, to study how galaxies fade and redden with time (the 
solid lines with crosses in Fig. 1). We construct a 'bright' 
(Mv-51og|,)/! --21) and 'faint' {Mv-5\og^Qh - -19) set of 
model galaxies. These models are not designed to be particu- 
larly realistic, rather they are discussed only to give a flavor for 
how star-forming galaxies may settle onto the red sequence if 
their star formation was stopped for some reason. We discuss 
the 'bright' set of model galaxies first, discussing the 'faint' set 
afterwards. The 'bright' set of model galaxies is based on a 
galaxy with a SFH of the form ■0 oc e"'' where is the star 
formation rate (SFR), t is the time elapsed since galaxy forma- 
tion, and T is the e-folding time of the SFH, set in this case to 
be 4 Gyr (this color would correspond roughly to a Sb galaxy 
in the local Universe). The formation redshift is assumed to be 
z/ = 2, and the metalUcity is assumed to be constant at solar, 
for simpUcity. The evolution of this model over cosmic time is 
shown by the bluest cross on the 'bright' track in each panel of 
Fig. 1. 

The sohd line at My - 5 logjQ /i ~ -21 in each panel of Fig. 1 
connects this r = 4 Gyr model with models which have the iden- 
tical SFH (crosses), except that they have had their SF stopped 
arbitrarily after an elapsed time of 90%, 80%, 50% and 10% of 
the age of the galaxy at the redshift interval of interest (these 
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Fig . 5 . — The luminosity density evolution of luminous red-sequence galax- 
ies. The luminosity density of those red-sequence galaxies destined to evolve 
passively into galaxies brighter than Mg — 51og|Q/j < -19.1 at (z) = 0.25 is 
shown by dashed, solid and dotted black lines for = 2. 3, 5 respectively. The 
smooth grey curves show the expected evolution of a galaxy population that 
completely formed at high redshift and simply aged to the present day. The 
error bars show uncertainties from cosmic variance. 

models correspond to the second bluest cross to the reddest 
cross, respectively; the bluest cross is the untruncated model 
as discussed above). When the SF is stopped, a burst of SF oc- 
curs which increases the stellar mass of the truncated model to 
the same stellar mass as the untruncated galaxy. The 'faint' set 
of models is analagous, except with 1/3 solar metallicity and 
a constant SFR (r = oo), in order to have the somewhat bluer 
colors characteristic of fainter galaxies. 

A cursory inspection of these tracks in Fig. 1 shows clearly 
that while faint red-sequence galaxies may plausibly be faded 
remnants of blue, later-type galaxies which exist at the same 
epoch, there are very few blue galaxies bright enough to fade 
into the most luminous red-sequence galaxies at all redshifts 
out to z ^ 1. Thus, luminous red-sequence galaxies cannot 
have originated from simple fading of single blue star-forming 
galaxies since redshift one. This implies that most luminous 
red-sequence galaxies were either completely formed earUer 
than z ~ 1 , or that these luminous red-sequence galaxies formed 
through mergers of galaxies. 

We explore these two possibilities by considering the lumi- 
nosity density in luminous red-sequence galaxies at z < 1 . Let 
us adopt as a null hypothesis that luminous red-peak galaxies 
form completely at high redshift and simply age to the present 
day (the first option). Then, if we choose to evaluate the lu- 
minosity density in red-sequence galaxies above a passively- 
fading luminosity cut (i.e., a constant stellar mass cut), we 
should observe a luminosity density which fades passively (like 
the grey passive evolution tracks in Figs. 4 and 5). If the total lu- 
minosity density in this sub-population does not follow a simple 
fading track, then a merger origin for at least some of the lumi- 
nous red galaxies at z < 1 appears necessary (the second option) 
to boost this population of luminous red galaxies. Recall that 
a 'bulge-i-fading disk' scenario is not allowed, as a reservoir of 
very luminous blue galaxies (blue because of their disks) would 



be required to fade and redden into red-sequence galaxies, in 
clear contradiction with the observations (Fig. 1). In order to 
place the most interesting constraints, we choose a luminosity 
cut corresponding roughly to the 'bright' My-51ogjo/j ^ -21 
track in Fig. 1 above which there are few blue galaxies which 
can fade onto the red sequence. 

The luminosity density of red-peak galaxies destined to pas- 
sively evolve into galaxies brighter than Mb - 5 logig h<-\9.l 
at (z) = 0.25 (as calculated using the rate of fading of Pegase 
model stellar populations in conjunction with an assumed for- 
mation redshift) is shown in Fig. 5 by dashed, solid and dotted 
black Unes for different formation redshifts z/ = 2 , 3 , 5 . This lu- 
minosity cut is equivalent to My-Slogjo/z < -20, accounting 
for the average rest-frame color of red-peak galaxies at that red- 
shift ofB-V = 0.9. We show example cuts as a function of red- 
shift in Fig. 3 as arrows, assuming Zf = 3. For reference, these 
example cuts areMs-Slogio/j < -19.1, -19.2, -19.32, -19.44, 
-19.54, -19.64, -19.74, -19.86, -20.01 for (z) =0.25, 0.35, 
0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.05, again for z/ = 3. We 
overplot large scale structure errors, as before. We show also 
the luminosity density in galaxies brighter than Mb— 5 logig h < 
-18.8 at z = for SDSS EDR red-peak galaxies as a sohd point; 
this cutoff corresponds to MB-5log^Qh < -19.1 at (z) = 0.25, 
passively faded to the present day. The passive fading expec- 
tation is shown in the grey lines. The luminosity densities are 
sums of the V/Vmax points, except for a tiny contribution to the 
two highest redshift bins, where we must extrapolate slightly 
using the Schechter function fit. 

It is clear that the observed luminosity density in luminous 
red peak galaxies evolves relatively Uttle over the interval < 
z < 1.1 (similarly to the luminosity density evolution in the 
whole red-sequence galaxy population). There is perhaps a 
modest increase out to z ~ 0.7, then a significant drop at higher 
redshift. This weak evolution does not agree with the expecta- 
tions of passive evolution models (where the galaxy population 
forms completely at high redshift and simply ages to the present 
day). Thus, the population of luminous red-peak galaxies (those 
destined to become red-peak galaxies with MB-51ogjg/z < 
-19.1 by (z) =0.25) grows in stellar mass by roughly a factor of 
two in the interval < z < 1.1. There is a hint that much of this 
stellar mass build-up happens at z > 0.7, with a modest gradual 
increase at later times. Clearly, however, a larger survey with 
lower uncertainties from large scale structure is required to fur- 
ther constrain the epoch of rapid stellar mass assembly in the 
luminous red galaxy population. 

The significance of the evolution in the luminous red galaxy 
population is that we have chosen a luminosity cutoff which 
corresponds to the 'bright' truncation model of Fig. 1. At all 
redshifts, it is clear that there are very few single blue galax- 
ies bright enough to fade into a red galaxy with luminosities 
brighter than this model. Yet, this luminous red galaxy popula- 
tion builds up in stellar mass by a factor of two in this redshift 
range. Red massive progenitors would have been included in 
the luminosity density estimates, and blue massive progenitors 
are observed to be very rare indeed (Fig. 1). Thus, merging 
between less luminous galaxies is required to build up these lu- 
minous red galaxies (see also, e.g., Le Fevre et al. 2000; Patton 
et al. 2002; Brinchmann & Ellis 2000, for discussion of the in- 
creased importance of merging in galaxy evolution at z ^ 1). 
Furthermore, owing to the rarity of luminous blue galaxies, 
these mergers must have at most a very brief luminous blue 
phase. This is possible if most mergers of massive galaxies are 
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either very dusty (obscuring and reddening the stellar popula- 
tions until the stellar populations age and redden sufficiently to 
join the red sequence) or are very gas poor to begin with (e.g., 
the elliptical-elliptical mergers observed by van Dokkum et al. 
2001, and explored theoretically by Khochfar & Burkert 2003). 
Exploring these two options in detail is far beyond the scope 
of this work, except to note that both have clear observational 
signatures. The highly-obscured gaseous mergers should be 
among the brightest far-infrared sources at any given redshift, 
and will therefore be well-constrained by current and planned 
observations in the far-infrared and sub- mm regime (e.g., Flores 
et al. 1999). On the other hand, the elUptical-elliptical mergers, 
if they build up the galaxy mass by more than a factor of two 
or so, will flatten the CMR at the highest luminosities (Bower, 
Kodama, & Terlevich 1998). Disentangling the contributions 
from both processes will be a fascinating prospect for the com- 
ing years. 

The distribution of galaxy colors also strongly constrains the 
'frosting' model of Trager et al. (2000). In this scenario, many 
early-types have residual low-level SF at z < 1. One could 
postulate that the factor of two evolution in stellar mass in 
red-sequence galaxies is due to residual SF in the other half 
of the galaxies, making them too blue to satisfy our redshift- 
dependent color cut. This residual SF slows to the present day, 
giving an apparent evolution of a factor of two in stellar mass 
while avoiding the need for wholesale galaxy assembly through 
galaxy mergers. However, this mechanism could not work for 
at least high-luminosity red galaxies. A 'frosting' of SF on top 
of a massive red-sequence galaxy at z ~ 1 would produce an 
even brighter blue galaxy, in clear contradiction with the obser- 
vations. Indeed, as estimated using the Pegase stellar popula- 
tion model, a color difference of A(U-V) ^ 0.25 (the color cut 
adopted in our definition of red-sequence galaxies) corresponds 
to an ongoing SF rate of only a few percent of the average SF 
rate at earUer times (i.e., flie birthrate b = SFR/(SFR) < 0.05). 
Thus, the observations strongly limit the amount of residual SF 
in massive red-sequence galaxies all the way from z ^ 1 to the 
present day to less than a few percent enhancement in stellar 
mass over the whole redshift range < z < 1.1. Our obser- 
vations say little about less luminous early-types; medium-to- 
high resolution, high S/N spectroscopy will allow us to more 
meaningfully constrain the SFHs of less luminous red-sequence 
galaxies. 

6.5. A Synthesis 

It is interesting to bring together some of the different threads 
in this paper into an overall qualitative, perhaps in parts naive, 
picture of galaxy evolution. 

Firstly, the color distribution is bimodal. As discussed above 
at the end of §6.4, birthrates (the relative present-day SFR com- 
pared to the average past SFR) for red-sequence galaxies are 
low < 5%. Interesting lower limits on the birthrates of blue- 
peak galaxies can be derived by assuming that these galaxies 
are dust-free. Again using the Pegase stellar population mod- 
els, the red end of the blue peak (with f/-y ~ 0.7) corresponds 
to ~ 0.1, whereas the bluest galaxies correspond to ~ 1 (us- 
ing Pegase models or the results of Kennicutt, Tamblyn, & 
Congdon 1994, to transform U-V into b) ^. Thus, galaxies 
on the blue peak need not be forming stars vigorously at the 
present day (see also Brinchmann et al. 2004). Indeed, modest 

^It should be noted that the 4th column of their Table 1 is mistakenly la- 
belled U-V, when it should be labelled U-B. 



rates of only 10% of the past average SFR are required to keep 
a galaxy on the blue peak (more, if the effects of dust are im- 
portant). Thus, bimodality in the color distribution of galaxies 
simply betrays the fact that young stars are very bright, and that 
if even a small number of young stars are present they domi- 
nate the optical colors of that galaxy. The red peak, in contrast, 
tells us that there is a substantial population of non-star forming 
galaxies, with present day SFRs of at most a few percent of their 
past average SFRs (for example, Hogg et al. 2002; Kauffmann 
et al. 2003; Bell et al. 2003, find that half or more of the stel- 
lar mass in the local Universe is in non-star forming early-type 
galaxies). 

Secondly, the B-band luminosity density of galaxies in the 
red sequence does not significantly evolve in the interval < 
Z < 1.1. This implies a build-up of stellar mass on the non-star- 
forming red peak by a factor of at least two since z ^ 1. This 
stellar mass must come from the blue peak galaxies, as we are 
not permitted to form the stars in situ because of the red colors 
of the red-sequence population. 

Is the relative number of galaxies in the red and blue peaks, 
and in the gap, consistent with this picture? We explore this 
by constructing a very simplistic model, in which some small 
fraction of star-forming galaxies in the blue peak are 'turned 
off (by some unspecified mechanism) and subsequently red- 
den and fade (see, e.g., Kodama & Bower 2001, for a model 
of this type appUed to rich galaxy clusters). Fading across the 
'gap' takes around ~ 40% of the Hubble time (as seen by the 
crosses on the fading tracks in Fig. 1). In order to reproduce 
the roughly factor of two to three difference in the number of 
galaxies in the gap compared to the blue peak, it is necessary 
to turn off roughly 5 - 10% of the blue galaxies per Gyr, letting 
them fade across the gap. This rate of transformation also gives 
a factor of roughly two increase in the stellar mass density on 
the red sequence. Thus, a very simpUstic picture where some 
small fraction of the blue population of galaxies have their SF 
stopped every Gyr (by some unspecified mechanism), subse- 
quently fading into red galaxies, holds water. 

Of course, we know that there must be other physical pro- 
cesses at work. A truncation-only model predicts that the char- 
acteristic luminosity of red galaxies is fainter than the blue 
population. Yet, the brightest and most massive present-day 
galaxies in the Universe are almost all red (e.g., Kochanek et al. 
2001; Bell et al. 2003). Furthermore, this situation persists out 
to z ~ 1 (Fig. 1). Another prediction is that the morphologies 
of the red and blue peak populations should be rather similar, 
which again is completely ruled out by observations (Strateva 
et al. 2001). Indeed, the hypothesis that luminous red-sequence 
galaxies at z = are faded remnants of bulge-l-blue disk systems 
at z ^ 1 whose disks have faded and reddened can be ruled out 
owing to the paucity of luminous blue galaxies at z ^ 1 . Galaxy 
mergers are a Ukely process for rectifying these shortcomings. 
Galaxy mergers will increase the characteristic luminosity of 
post-merger galaxies through the addition of the pre-existing 
stellar populations, in conjunction with any SF associated with 
the merger itself. Furthermore, the violent relaxation associated 
with galaxy mergers leads to more spheroidal-dominated mor- 
phologies (e.g., Eggen, Lynden-BeU, & Sandage 1962; Toomre 
& Toomre 1972; Barnes & Hemquist 1996). Indeed, mergers 
and galaxy interactions may even provide a mechanism for the 
truncation of SF, through the consumption of the gas in an in- 
duced starburst (e.g., Barnes & Hemquist 1996; Barton, Geller, 
& Kenyon 2000). The only constraint which we can place on 
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the role and frequency of these mergers with these data is that 
the mergers responsible for the most luminous red galaxies can- 
not have a long blue phase, owing perhaps to dust or a lack of 
gaseous content in the progenitors. 

This qualitative picture agrees well with many of the features 
of current models of galaxy formation in a hierarchical CDM 
Universe (e.g., White & Frenk 1991; Kauffmann et al. 1996; 
Somerville & Primack 1999; Cole et al. 2000). Furthermore, 
it makes a number of predictions that are testable in the near 
future. For example, the merger origin of the most luminous 
red galaxies can be tested by both deep optical HST imaging 
over wide fields and by far-infrared and sub-mm deep imaging 
to constrain the frequency of mergers between gas-poor early 
types and more gaseous later types, respectively (building on 
the studies of, e.g., Flores et al. 1999; Le Fevre et al. 2000). 
Additionally, one can test the viability of the truncation model, 
and constrain truncation mechanisms, through a combination of 
deep HST imaging (for galaxy morphologies as a function of 
color and magnitude) and deep spectroscopy (to check for sig- 
natures of recently- truncated SF, such as enhanced B aimer ab- 
sorption lines), extending the cluster-based work of, e.g., Pog- 
gianti et al. (1999) to the more general field environment. Fur- 
thermore, one can seek to explore these different galaxy popu- 
lations in different local environments, to seek the signature of 
more rapid environmentally-induced transformation in denser 
environments, extending to higher redshift the important work 
of, e.g., Lewis et al. (2002), Gomez et al. (2003), and Balogh et 
al. (2004). 

7. CONCLUSIONS 

The evolution of red- sequence/early-type galaxies directly 
reflects the importance and number of major galaxy mergers, 
and therefore strongly constrains hierarchical models of galaxy 
formation and evolution. We have explored the rest-frame col- 
ors and luminosities of 25000 R <24 galaxies from 0.78 
square degrees from the COMBO- 17 survey, encompassing a 
cosmologically-representative total volume of ^ lO^ft-^Mpc^. 
We find that the distribution of galaxy colors is bimodal at all 
redshifts out to z ~ 1. 

The blue star-forming peak has colors that become redder to- 
wards the present day at a given magnitude, indicating changes 
in characteristic stellar age, metallicity or dust content with 
epoch. Furthermore, there are many more luminous blue galax- 
ies at z > 0.5 than there are at the present day (this is discussed 
in much more detail by Wolf et al. 2003). 

The red non-star-forming peak forms a scattered but well- 
defined color-magnitude relation (CMR) at0.2<z< 1-1, with 
luminous galaxies being redder than their fainter counterparts. 
The color of red galaxies at a given rest-frame magnitude be- 
comes bluer with increasing redshift. The quantitative size of 
this change is consistent with passive aging of ancient stellar 
populations to the present day. 

Using this empirically-motivated definition of early-type gal- 
axies, we estimate the rest-frame fi-band luminosity function 
and luminosity density of red-sequence galaxies in the inter- 
val < z < 1.1. There may be some contribution from dusty 
star-forming galaxies at the highest redshifts, making our re- 
sults upper limits. Nevertheless, we find mild evolution of the 
rest-frame B-band luminosity density between < z < 1.1. An- 
cient stellar populations would fade by a factor of two to three 
in this time interval; therefore, this mild evolution betrays an 
increase in the steUar mass on the red sequence since z ~ 1 by 
at least a factor of two. This evolution is consistent qualitatively 



and quantitatively with the evolution expected from the hierar- 
chical build-up of stellar mass via galaxy mergers in a ACDM 
Universe. The largest source of error is large-scale structure, 
implying that considerably larger surveys are necessary to fur- 
ther refine this result. 

Finally, we explore the evolution of the red, blue and gap 
population on the color-magnitude plane. We find that a sce- 
nario in which SF is stopped (for an as-yet-unspecified reason) 
in 5-10% of the blue galaxy population per Gyr reproduces 
the relative numbers of red, blue and gap galaxies. We argue 
also, based on the luminosity functions and morphologies of 
red galaxies in the local and distant Universe, that galaxy merg- 
ing plays an important role. In particular, at least some of the 
the most luminous red galaxies must be formed in galaxy merg- 
ers which are either very dusty or take place between gas-poor 
progenitors. Disentanghng further the processes which drive 
the evolution of the red galaxy population at z < 1 wiU indeed 
present a fascinating observational chaUenge over the next few 
years. 
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APPENDIX 



RED-SEQUENCE GALAXIES IN THE LOCAL UNIVERSE 

Because of its homogeneity and color information, the SDSS 
Early Data Release (EDR; Stoughton et al. 2002) offers an ideal 
comparison sample to COMBO- 17 at redshift zero. In order to 
compare as fairly as possible with COMBO- 17, we choose to 
synthesize U, B, and V magnitudes from SDSS ugr data, using 
the UBV bandpasses that were used to construct COMBO-17 
rest-frame magnitudes. We adopt the transformations: 

f/-y =-0.713 -I- 0.826(M-r);cr- 0.023, (Al) 

g-B = -0.155-0.370(^-r);a~ 0.012, (A2) 

^-y = -0.010-l-0.609(g-r);CT~ 0.006, (A3) 

which were derived from the single stellar populations from the 
Pegase stellar population synthesis model, and are vaUd for 
stellar populations between 1/200 and 5 x solar metallicity, and 
ages between 50Myr and 20Gyr These transformations are 
dominated by the systematic zero point uncertainties (in ab- 
solute terms) of the SDSS and Johnson photometric system, 
and by limitations of the SDSS Petrosian magnitude. These 
uncertainties are likely < 0.1 mag in terms of color, although 
absolute calibration of u or t/-band data is always challeng- 
ing. SDSS EDR galaxies were chosen to have measured red- 
shifts and a reddening-corrected Petrosian magnitude in r be- 
tween 14.5 < r < 16.5, yielding a total sample of 3592 galaxies. 
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Fig. A6. — Rest-frame synthesized U-V colors and B-band luminosity function of SDSS EDR galaxies. The left-hand panel shows the synthesized U -V vs,. My 
color-magnitude diagram for a volume-weighted sample of 14.5 < r < 16.5 SDSS galaxies. We add random offsets of (J(ij-v) = 0.01 and ay =0.1 for clarity, thus 
giving clumps of data points for galaxies which are included multiple times owing to their small accessible volume in SDSS. The solid and dotted lines show the 
best-fit CMR and the Butcher-Oemler-style cutoff adopted in this paper, assuming a slope of -0.08 (following Bower, Lucey, & EUis 1992; Terlevich, Caldwell, & 
Bower 2001). The inset panel shows the color distribution of these galaxies, with the slope of the CMR taken out. The right-hand panel shows the rest-frame B-band 
luminosity fijnction of the red-sequence galaxies, assuming a 0.3 mag magnitude error caused by photometric redshift uncertainty. The solid circles with error bars 
show the V/Vmax estimates, and the solid Une a Schechter fit to these estimates. The dotted line shows the Schechter fit to the B-band luminosity function without 
the additional 0.3 mag scatter. 



The Petrosian magnitude of concentrated Cr > 2.6 (i.e. early- 
type) galaxies is corrected by -0.10 mag to crudely account for 
the well-known shortfall of Petrosian magnitudes for centrally- 
concentrated galaxies (e.g., Blanton et al. 2003b). 

We show the distribution of synthesized U-V color against 
My for a volume-limited SDSS subsample in the left-hand 
panel of Fig. A6. Galaxies are weighted by the inverse of the 
volume within which they are observable (1 /Vmax), and are se- 
lected to have a uniform distribution in terms of 1 / Vmax- Since 
SDSS is an apparent magnitude-limited survey, only a small 
fraction of the bright, distant galaxies are selected for inclusion 
on this plot, whereas the fainter galaxies which are visible over 
a much smaller volume are selected multiple times. We choose 
to construct a sample of 2000 galaxies with My - 5 logjQ /i < - 1 8 
(to roughly match the magnitude cutoff for galaxies at z ^ 0.6), 
and we add random color and magnitude offsets with a(u-v) = 
0.01 and civ = 0.1 to allow the reader to see multiply-included 
galaxies explicitly. This SDSS control sample shows a CMR 
with scatter 0.16 mag; this large scatter is expected, as the to- 
tal u-r colors we used have typical uncertainties ^0.15 mag, 
after observational error and evolution and TT-correction errors 
are accounted for. The intercept of the CMR is consistent with 
the intercept of the Abell 754 galaxy cluster (Mcintosh, Rix, 
& Caldwell 2004) and the Nearby Field Galaxy Survey (Jansen 
et al. 2000) morphologically-classified early-type galaxy CMR 
to within 0.1 mag. This kind of offset is not unexpected given 
the difficulty of caUbrating absolute photometry to better than 
10%. 

Following Bell et al. (2003), we construct a rest-frame B- 
band luminosity function of red-sequence galaxies (redder than 
the dashed line in Fig. A6) using the V /Vmax formalism. Galaxy 
selection for inclusion in this luminosity function is unbiased by 
any zero-point uncertainty in caUbration, as this criterion is ap- 
plied relative to the ridge-line of the CMR, which is measured 
using the same data. In order to constrain the effects of magni- 
tude error (from photometric redshift error) on the luminosity 
function (as discussed by Chen et al. 2003), we convolve our B- 



band luminosities with a random 0.3 mag error, which simulates 
a 5z ^ 0.02 at z ^ 0.25 (our lowest redshift bin). We also con- 
strain the faint-end slope a = -0.6, to match with the COMBO- 
17 luminosity function determinations. The Schechter func- 
tion fit parameters are cj)* = 10.2 ±0.9 x 10"^/i^Mpc~^mag"\ 
M* - 51og|o h = -19.16 ± 0.10, a = -0.6, and Jb = 6.4 ± 0.2 x 
IO^Lq/iMpc"^ (random errors only). The luminosity function 
calculated in this way is shown by the sohd black circles and 
solid line in the right-hand panel of Fig. A6. For reference, a 
Schechter function fit to the unconvolved (i.e. real) B-band syn- 
thesized luminosity function is shown by the dotted line. The 
M* of this fit is ~ 0.05 mag fainter than the convolved fit (with 
corresponding increase of (j)* to compensate); this modest offset 
is the maximum bias that we expect in the Schechter fit param- 
eters because of our use of 17-passband derived photometric 
redshifts. This is in agreement with Chen et al. (2003), who 
show that we would expect only a modest bias towards brighter 
M* because of our relatively small photometric redshift uncer- 
tainty. The uncertainty in luminosity density is roughly 10% in 
a systematic sense; this error is dominated by uncertainties in 
absolute magnitude calibration and in correcting the luminosity 
function normalization to the observed all-sky \b\ > 30° density 
of 10</:< 13.5 galaxies. 

For more discussion of the selection, photometric parame- 
ters, V /Vmax luminosity function estimation, and uncertainties, 
see Bell et al. (2003). 



